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Abstract. We present broad-band V and / imaging and long-slit spectroscopy in the optical range AA3600 - 7500A 
of the dwarf irregular galaxy HS 1442+4250. The oxygen abundance 12 + log(0/H) = 7.63 ± 0.02 (Z = Z e /19)* 
in the brightest H II region of HS 1442+4250 places the galaxy among the most metal-deficient emission-line 
galaxies. The low metallicity and blue colour (V - I) ~ 0.4 mag of the low-surface-brightness (LSB) component 
make HS 1442+4250 a likely rare young dwarf galaxy candidate. We use four methods to estimate the stellar 
population age in the LSB component of HS 1442+4250. Different star formation histories are considered. The 
equivalent widths of hydrogen Ha and H/3 emission lines, and of hydrogen H7 and H<5 absorption lines, the spectral 
energy distribution and the observed (V — I) colours of the LSB regions are reproduced quite well by models with 
only young and intermediate-age stellar populations. By contrast, the observational data cannot be reproduced 
by a stellar population formed continuously with a constant star formation rate in the age range from to <; 
2 Gyr. While a faint old stellar population in HS 1442+4250 with an age ii 2 Gyr is not excluded, we find no 
evidence for such a population from the present data. 
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1. Introduction 

The dwarf irregular galaxy HS 1442+4250 (= UGC 9497) 
from the Hamburg Survey (HS) was first classified as an 
emission-line galaxy by Sanduleak & Pesch l|1982fl . With 
the coordinates a(J2000.0) = 14 h 44 m 12?l, <5(J2000.0) = 
+42°37'37" it is situated in the direction of a low-density 
region in the galaxy spatial distribution. HS 1442+4250 
has been studied by Tifft et al. I|1986[) . Popescu et al. 
(1990|, and Popescu & Hopp (2000) who were searching 
for dwarf galaxies in voids. Popescu & Hopp H2000fl have 
measured the fluxes and equivalent widths of the emis- 
sion lines in HS 1442+4250. They derived its redshift as 
z = 0.0025 and an oxygen abundance 12 + log(0/H) = 
7.89. The galaxy has also been studied spectroscopically 
by Kniazev et al. (|1998[) and Pustilnik et al. I|1999[) . who 
derived 12 + log(0/H) = 7.68 and 7.7, respectively. 

B and R surface photometry of HS 1442+4250 has 
been presented in Vennik et al. I|1996fl . This galaxy with 
Mb = -15.2 mag was found to be one of the bluest dwarf 



galaxies in their sample with (B — R)^ = 0.42 mag. This 
implies that the light from HS 1442+4250 is dominated 
by a relatively young stellar population. 



Send offprint requests to: N.G. Guseva, guseva@mao.kiev.ua 
* 12 + log(O/H) = 8.92 (Anders & Grevesse \TM)l . 



In this paper we present V and I photometric and 
spectroscopic observations of HS 1442+4250. We study 
its properties and put constraints on the age of its low- 
surface-brightness (LSB) component. Throughout this pa- 
per we adopt a distance of 12.4 Mpc for the dwarf galaxy 
derived from its redshift z — 0.0025 and a Hubble constant 
of 75 km s _1 Mpc -1 , after correction for Virgocentric 
infall (Kraan-Korteweg [1986). At this distance 1" cor- 
responds to a linear size of 60 pc. The structure of the 
paper is as follows. In Sect. El we describe the observa- 
tions and data reduction. The photometric properties of 
HS 1442+4250 arc described in Sect. EJ In Sect. H we de- 
rive chemical abundances in its two brightest H 11 regions. 
The properties of its stellar populations are discussed in 
Sect. |51 Finally, Sect. summarises our main conclusions. 
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Fig. 1. V image of HS 1442+4250. The overlayed V contours correspond to surface brightness between 20.5 mag 
arcsec -2 and 26.5 mag arcsec -2 in steps of 0.5 mag. 



2. Observations and data reduction 

2.1. Photometric observations 

Broad-band V and I images were obtained with the Kitt 
Peak 2.1m telescope 1 on April 18, 1999 during photo- 
metric conditions. The telescope was equipped with a 
thinned Tektronix 1024 x 1024 CCD operating at a gain 
of 3e - ADU -1 and giving an instrumental scale of 0'.'305 
pixel -1 and a field of view of 5'. The total exposures of 
20 min in V and 40 min in / were split into three subex- 
posures, slightly offset with respect to each other for re- 
moval of cosmic particle hits and bad pixels. The point 
spread function in V and / were respectively 1'.'78 and 
1"67 FWHM. 

Bias- and flat-field exposures were taken during the 
same night. The data reduction was done using IRAF 2 . 
Images were calibrated by observing four different stan- 
dard fields from Landolt (1992), each 3-4 times at differ- 
ent airmasses during the same night. Our calibration un- 

1 Kitt Peak National Observatory (KPNO) is operated by 
the Association of Universities for Research in Astronomy 
(AURA), Inc., under cooperative agreement with the National 
Science Foundation (NSF). 

IRAF is the Image Reduction and Analysis Facility dis- 
tributed by the National Optical Astronomy Observatory, 
which is operated by the AURA under cooperative agreement 
with the NSF. 



certainties are estimated to be well below 0.05 mag in all 
bands. Reduction steps included bias subtraction, removal 
of cosmic particle hits, flat-field correction and absolute 
flux calibration. 

2.2. Spectroscopic observations and data reduction 

The spectroscopic observations of HS 1442+4250 were car- 
ried out on June 18, 1999, with the Kitt Peak 4m Mayall 
telescope, in combination with the Ritchey-Chretien spec- 
trograph and the T2KB 2048 x 2048 CCD detector. The 
slit was centered on the brightest star-forming region and 
oriented along the elongated body of the galaxy with posi- 
tion angle PA. = 60?4, close to the direction of its major 
axis (see Fig.^l. The slit orientation during the observa- 
tions was close to the parallactic angle to minimize the 
effects of differential refraction. Hence, no correction was 
made for this effect. A 2" x 300" slit with the KPC-10A 
grating in first order and a GG 375 order separation fil- 
ter were used. The spatial scale along the slit was 0'.'69 
pixel -1 and the spectral resolution ~7 A (FWHM). The 
spectra were obtained at an airmass of 1.44. The total ex- 
posure of 60 minutes was broken up into 3 subexposures. 
Two Kitt Peak spectrophotometric standard stars, Feige 
34 and HZ 44, were observed for flux calibration. Spectra 
of a He-Ne-Ar comparison lamp were obtained for wave- 
length calibration. 
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Fig. 2. (a) Contrast-enhanced I contour map of HS 1442+4250. The brighter regions labeled c and e and the fainter 
regions a, b, d, /, g and h are arranged along the major axis of the galaxy, over a projected length of ~ V (3.5 kpc). 
Triangles, labeled I and II, show the positions of the outermost regions with Ha and H/3 emission and without any 
absorption features. Sources labeled 1 to 5 are significantly redder (0.8 < (V — I) ^ 2.0) than the LSB component 
((V — I) ~ 0.4 mag) and most probably are not related to HS 1442+4250. (b) V surface brightness distribution along 
the major axis of HS 1442+4250 computed within a rectangular strip of 2" width. The origin is at the location of 
region c. The six compact regions indicated in the upper panel from a to / span a range in surface brightness between 
20 and 22.6 V mag arcsec -2 and regions g and h between 22.3 and 23 V mag arcsec -2 . (c) (V — I) colour distribution 
within a 2" strip along the major axis of HS 1442+4250. The colour between regions a and h is nearly constant at 
~ 0.4 mag (dotted line). The bluest region ((V — I) ~ -0.5 mag) is coincident with region c. Region a presents a 
large spatial colour variation and contains a compact source which is much fainter on the V image. This source has 
a maximal (V — I) colour of 0.69 mag, being significantly redder than the adjacent regions (see also Fig. 121. Filled 
circles show the colours predicted for a stellar population formed continuously with a constant star formation rate in 
two periods with 1) intermediate-age and 2) young stellar population. The filled circle for region c indicates the colour 
of a 4 Myr stellar population to which has been added the observed ionized gas emission in region c. The asterisk 
shows the colour of a 4 Myr pure stellar population emission, while the open triangle shows the colour of ionized gas 
emission. 



The data reduction was made with the IRAF software 
package. This includes bias subtraction, flat-field correc- 
tion, cosmic-ray removal, wavelength calibration, night 
sky subtraction, correction for atmospheric extinction and 
absolute flux calibration of the two-dimensional spectra. 

One-dimensional spectra for abundance determination 
in the two brightest H n regions c and e (Fig. EJ) were 
extracted within apertures of 2" x 2'.'8. 

In addition we extracted one-dimensional spectra of 
six regions in the LSB component of the galaxy showing 



hydrogen Balmer absorption lines. The selected regions, 
labeled in Fig. a, b, d, f, g and h, are listed in Tables^] 
and [3 together with their positions relative to the bright- 
est region c and their angular extent along the slit. We 
also extracted one-dimensional spectra of two additional 
outermost regions labeled I and II with only Ha and H/3 
emission. 
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3. Imaging 

3.1. Morphology 

Photometric signatures of a young to moderately evolved 
stellar population are present all along the major axis of 
HS 1442+4250. The (V - I) colour map (Fig. EJ shows 
that the bluest region coincides with the brightest region 
c (Fig.EJ where (V - I) ~ -0.5 mag. This region is partly 
resolved and has an effective radius of 2'.'4 («140 pc) and 
a mean FWHM of 2'.'5 in V. Its apparent V magnitude 
of 18.0 mag corresponds to an absolute V magnitude of 
~ -12.5 mag. The (V — I) colour of the second brightest 
region e (V = 19.65 mag) is ~ 0.2 mag, markedly bluer 
than the local LSB component with an average (V — I) 
colour of ~ 0.4 mag. Sources c and e contribute more than 
80% of the flux in excess of the LSB component inside the 
25 V mag arcsec -2 isophote. Several other fainter compact 
sources are seen along the major axis of the galaxy (Fig. 
Et)- Variations of the (V — I) colour on spatial scales of 
2" - 3" (Fig. [3J imply that non-uniform extinction may 
be significant. Alternatively, it may be that part of the 
reddest small regions are due to background/foreground 
red sources. This is likely to be the case for the compact 
region a, clearly detected in both the / image and the 
colour map, but marginally seen in V. Roughly 15" SE of 
region c there is a red, relatively bright (V = 18.6 mag) 
source labeled 5 with a mean (V — I) i; 1.2 mag, i.e., 
roughly 0.6 mag redder than the LSB component of HS 
1442+4250 (see Fig.0jD). The surface brightness profile of 
this source can be fitted by a Sersic profile (Sersic EEEJ 
with a V central surface brightness of rj 21.9 mag arcsec -2 
and an exponent 77 of ~ 0.7 (see Eq. 1 in Guseva et al. 
2003 ). Most likely, this source as well as those labeled 1 - 

4 in Fig. 01 are background galaxies. 

3.2. Surface brightness and colour distribution 

The surface brightness profiles (SBPs) of HS 1442+4250 
in Fig. EH were obtained using one of the methods (method 
iii) described in Papaderos et al. I)1996a|l and Guseva et al. 
(20031 after the red compact sources 1 through 5 in Fig. 
|2H have been fitted and subtracted from the combined 

V and I images. The extended sources 4 and 5 overlap 
with the LSB component of HS 1442+4250 at a surface 
brightness fainter than 24.5 V mag arcsec -2 . In this case, 
a two-dimensional model was fitted to their light distribu- 
tion and subtracted from the original images. To ensure 
that possible residuals in the subtraction of sources 4 and 

5 do not affect the surface photometry, SBPs were com- 
puted by fitting ellipses to isophotes after screening out 
the area subtended by position angles -145° < < -95° 
with respect to region c. 

The SBPs of HS 1442+4250 (Fig. Hi) are well fitted 
by an exponential in the radius range 15"^ R* Ss 26", 
corresponding to a surface brightness fainter than ~ 24 

V mag arcsec -2 . In this range we obtain a V band scale 



length of w 4" 5, in good agreement with the value of 4" 2 
derived previously in the B band by Vennik et al. J2000). 

However, an inward extrapolation of the exponential 
fit predicts for radii between 5" and 10" a higher intensity 
than that observed (not shown in Fig. HJi)- Therefore an 
exponential model flattening for small radii appears to 
be necessary to adequately fit the LSB emission of HS 
1442+4250. This conclusion is in line with the results by 
Vennik et al. (2000) who noted the convex shape of the 
SBP of HS 1442+4250 over a substantial radius range. 
Indeed, for radii > 15", by fitting a Sersic ljl968|l profile 
to the SBPs in Fig. 0] we obtain an exponent 77 of ~ 1.2 
and 1.3 in V and /, respectively, slightly greater than that 
corresponding to a pure exponential distribution (77 = 1). 

Alternatively, a fitting formula which reproduces an 
inward flattening of the exponential distribution has been 
discussed by Papaderos et al. I|1996al their Eq. 22; see 
also Guseva et al. 120031 for a detailed explanation) . Near 
the center, such an intensity distribution depends on the 
relative central intensity depression q = AI/Iq, where Iq 
is the central intensity obtained by extrapolation of the 
outer exponential slope to R* = 0, and on the cutoff- 
radius ba inside of which the central flattening occurs. 

We find that fits to the SBPs for radii R* > 15" adopt- 
ing a Sersic law are nearly indistinguishable from those 
given by Eq. 22 of Papaderos et al. I|1996af) with b — 2.1 
and q = 0.3 (Fig. Both fits predict a central surface 
brightness of ~ 21.8 V mag arcsec -2 for the LSB compo- 
nent (Fig. 2H)j which is roughly 1.3 mag fainter than the 
value fj,E,o predicted by extrapolation of the exponential fit 
from the outer SBP part to R* = 0. The surface brightness 
distribution of the residual emission in excess of the LSB 
model (Fig.^Hj crosses) is to be attributed to star- forming 
regions along the major axis of HS 1442+4250 (starburst 
component). The integrated V magnitudes inside the 25 
and 27 V mag arcsec -2 isophotes are 17.65 and 17.35 mag, 
respectively. Such values imply that the starburst compo- 
nent contributes & 20% and 10% of the HS 1442+4250 
total light in the V inside the respective isophotal radii. 
About one half of the total starburst component emission 
originates from region c (Fig. |2H, EH) ■ 

Table 1 summarizes the derived photometric quanti- 
ties. They are not corrected for inclination and foreground 
Galactic extinction. Cols. 2 and 3 give, respectively, the 
central surface brightness /ze,o and scale length a of the 
LSB component obtained from exponential fits to the 
SBPs for R* > 15" and weighted by the photometric un- 
certainty of each point. For the apparent major-to-minor 
axis ratio of 4.1 at 26.5 V mag arcsec -2 and assuming 
a disk geometry for the LSB component the inclination- 
corrected central surface brightness is ~ 1.2 mag fainter 
than the value quoted in col. 2. In columns 4 and 5 we 
show the apparent magnitudes ttilsb of a pure exponen- 
tial distribution (q = 0) and m^ SB of an inwards flatten- 
ing distribution (Eq.22 in Papaderos et al. 1996a) with the 
derived flattening parameters (b,q) = (2.1,0.3). Note that 
both distributions refer to the same extrapolated central 
surface brightness he,o and a (cols. 2 and 3). Columns 6 
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Fig. 3. (V — I) colour map of HS 1442+4250 displayed in the range 0.25 to 0.55 mag. The colours are corrected for 
foreground Galactic extinction (Ay = 0.044 mag, Aj — 0.026 mag). The overlayed contours are from 20.5 to 26 V 
mag arcsec -2 in steps of 0.5 mag. Crosses along the major axis of the galaxy mark the positions of regions a, / (black 
crosses) and c, e (white crosses) (see Fig. EK)- The positions of the red sources 4 and 5 are marked by squares. The 
bluest region with (V — I) ~ -0.5 mag coincides with region c. 



Table 1. Structural properties of the starburst and LSB components. 



Band 


A*E,0 


a 




m LSB 


P 25 


mp 25 


E25 


m E25 


msBP 


TOtot 


r e s,rso 




mag arcsec -2 


pc 


mag 


mag 


pc 


mag 


pc 


mag 


mag 


mag 


pc 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


V 


20.50+0.08 


272+ 6 


15.23 


15.54 


496 


17.65 


1127 


15. 65 


15.37+0.02 


15.35 


503,822 


I 


20.18+0.20 


291+15 


14.76 


15.07 


709 


17.31 


1290 


15.16 


14.95+0.03 


14.93 


556,880 



through 9 list quantities obtained from profile decompo- 
sition whereby the intensity distribution of the LSB com- 
ponent was modeled by the modified exponential distribu- 
tion proposed by Papaderos et al. (1996a). In cols. 6 and 
8 we tabulate the radial extents P25 and E25 of the star- 
burst and LSB components respectively, both determined 
at a surface brightness 25 mag arcsec -2 (see Fig.0jD). The 
respective apparent magnitudes of each component within 
P25 and E25 are listed in cols. 7 and 9. The total magnitude 
of HS 1442+4250 derived from SBP integration out to the 
outermost radius of « 26" is listed in col. 10. The appar- 
ent magnitude of the galaxy, obtained within a polygonal 
aperture after removal of sources 1 through 5 (Fig. Et) 
is given in col. 11. In col. 12 we list the effective radius 
r c g and the radius r 80 which encircles 80% of the galaxy's 
total flux. 

As shown in Fig. the (V — I) colour of HS 
1442+4250 beyond the radius P25 increases with an av- 
erage gradient of 0.24 mag kpc -1 . The observed gradual 
colour increase is in agreement with the linear slope de- 
rived by subtracting the modeled distributions for the LSB 
component (thick-gray lines in Fig.^J}). From the colour 
profile we derive a mean (V — I) = 0.63 ± 0.05 mag in the 



radius range 15" < R* < 26". The (V - 1) colour at R* > 
18" is uncertain, and it is not clear whether the colour gra- 
dient is attributable to HS 1442+4250, or caused by the 
background sources and residuals after their subtraction. 

4. Chemical abundances 

In this Section we analyze the element abundances in the 
two brightest H 11 regions c and e. Their spectra, shown in 
Fig. are characterised by strong nebular emission lines. 

The emission line fluxes were measured by Gaussian 
profile fitting. The errors of the line fluxes include the er- 
rors in the fitting of profiles and those in the placement 
of the continuum. They have been propagated in the cal- 
culations of the elemental abundance errors. The spectra 
were corrected for interstellar extinction with the extinc- 
tion coefficient C(H/3) derived from the hydrogen Balmer 
decrement using the equations from Izotov et al. ([1994) 
and the theoretical hydrogen emission line flux ratios from 
Brocklehurst l|1971[) . The emission line fluxes relative to 
the H/3 emission line flux, both observed (F(A)) and cor- 
rected for extinction and underlying stellar hydrogen ab- 
sorption (/(A)), the equivalent widths EW of the emission 
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lines, the observed fluxes of the H/3 emission line, and the 
equivalent widths of the hydrogen absorption lines for re- 
gions c and e are listed in Tabled 

Our measured fluxes of the emission lines in region c 
are in fair agreement with those from Popescu & Hopp 
(2000) except for the [O n]A3727A emission line which in 
our case is ~ 2 times weaker. The agreement between our 
measured emission line fluxes and those by Kniazev et al. 
lfT5MI> and Pustilnik et al. iTlM5| is much better. 

Note that the flux of the He II A4686 nebular emission 
line in region c is relatively high (~ 3% of the H/3 flux) 
(Table EJ implying that hard radiation beyond A228A is 
intense. Such a property is not unusual for hot supergiant 
H II regions in low-metallicity dwarf galaxies. Strong He II 
A4686 nebular emission has also been detected in I Zw 18 
(Z Q /50, e.g ., Izotov et al. lT555l> . SBS 0335-052 (Z /4O, 
Izoto vet al.[TW7\\mm . Tol 1214-277 (Z Q /23, e.g., Fricke 
et al. HSU), SBS 1415+437 (Z Q /21, Thuan et al. 175351 
1999J and some other low-metallicity galaxies (Guseva et 

ai. rennnii . 

Only a few strongest emission lines are detected in the 
spectra of the two outermost regions I and II. The ob- 
served fluxes and EWs of the Ha and H/3 emission lines 
in these regions are listed in Table Since the continuum 
in the red part of the spectra in those regions is weak, the 
EW(Ra) are very uncertain. The extinction coefficient 
C(H/3) for regions I and II is set to 0. 

The photoionization model used to convert line fluxes 
into abundances is the same as described in Guseva et al. 
(2003}. The ionic and heavy element abundances for re- 
gions c and e together with electron temperatures, electron 
number densities and adopted ICFs are given in Table 

4.1. Heavy element abundances 

The oxygen abundances in regions c and e are 12 + 
log(0/H) = 7.63 ± 0.02 (Z ~ Z Q /19) and 7.54 ± 0.20 {Z 
~ Zq/24), respectively. These two determinations agree 
within the errors. However, the oxygen abundance in re- 
gion e may have been underestimated. This is because the 
weak [O m]A4363 emission line [7(A4363)/J(H/3)=0.055 
(TableEJ)] may be enhanced by shocks. This effect is likely 
lower in region c where the [O lii]A4363 emission line is 
stronger. For comparison, Popescu & Hopp (2000) ob- 
tained for HS 1442+ 4250 12 + log(0/H) = 7.8 9, while 
Kniazev et al. I|1998JI and Pustilnik et al. (1999) derived 
12 + log(0/H) = 7.68 and 7.7, respectively. The former 
value is significally larger than ours due to the higher 
[O n]A3727A and lower [O m]A4363A fluxes measured by 
Popescu & Hopp PTOjl . 

The neon-to-oxygen abundance ratio log Ne/O = -0.75 
for region c is in good agreement with the mean ratio de- 
rived from previous studies of BCDs (e.g., Izotov & Thuan 
1999f). The nitrogen-to-oxygen abundance ratio log N/O 
= -1.44 is ~ 0.15 dex higher than the N/O ratios obtained 
by Thuan et al. ifTMB) and Izotov & Thuan HMty for the 
very metal-deficient BCDs with Z ^ Z Q /20. 



Table 3. Element abundances in regions c and e. 



Value 


region c 


region e 




17810+320 


16410+4100 


T c (0 II) (K) 


14960+250 


14390+3420 


T C (S IIl)(K) 


16480+260 


15320+3400 


iV c (S ll)(cm- 3 ) 


10+10 


10+10 


O+/H+(xl0 5 ) 


0.467+0.022 


1.372+0.890 


O +2 /H+(xl0 5 ) 


3.685+0.159 


2.067+1.279 


O +3 /H+(xl0 6 ) 


1.464+0.194 




O/Hf' x\0 5 ) 


4 299+0 188 


3 439+1 558 


12 + log(0/H) 


7.63+0.02 


7.54+0.20 


N+/H+(xl0 7 ) 


1.702+0.128 


5.361+2.585 


ICF(N) a 


9.20 


2.51 


lUg^l\/ \J J 


-1 438+0 038 


-1 408+0 287 


Ne +2 /H + (xl0 5 ) 


0.662+0.035 


0.467+0.300 


ICF(Ne) a 


1.17 


1.66 


log(Ne/0) 


-0.746+0.030 


-0.646+0.342 


S+/H+(xl0 7 ) 


0.700+0.036 




S +2 /H+(xl0 7 ) 


5.593+0.725 




ICF(S) a 


2.30 




log(S/0) 


-1.473+0.054 




Ar+ 2 /H+(xl0 7 ) 


1.185+0.070 




Ar+ 3 /H+(xl0 7 ) 


1.420+0.252 




ICF(Ar) a 


1.01 




log(Ar/0) 


-2.213+0.048 




He+/H+(A4471) 


0.075+0.007 




He+/H+(A5876) 


0.078+0.003 




He+/H+(A6678) 


0.077+0.006 




He+/H+ (mean) 


0.078+0.002 




He +2 /H+(A4686) 


0.003+0.000 




He/H (A4471) 


0.078+0.007 




He/H (A5876) 


0.081+0.003 




He/H (A6678) 


0.080+0.006 




He/H (mean) 


0.080+0.002 




Y (A4471) 


0.237+0.022 




Y (A5876) 


0.244+0.009 




Y (A6678) 


0.242+0.018 




Y (mean) 


0.243+0.008 





^ICF is the ionization correction factor. 



4.2. He abundances 

The He I emission lines in region c are strong, allowing 
a reliable determination of the He abundance. Five He I 
A3889, A4471, A5876, A6678, A7065 emission lines are used 
to correct their fluxes self-consistently for collisional and 
fluorescent enhancement (see the method description in 
HI 



Izotov et al. 199 IJ. Then the He + abundance is derived 
from the corrected He I A4471, A5876, A6678 emission line 
fluxes. Some fraction of He (~ 3.8%) is in the Hc +2 form 
as derived from the He n A4686 emission line flux. The 
total helium abundance is derived as the sum of the He + 
and He +2 abundances and is shown in Table [3] for each of 
the three He I lines used for the abundance determination. 
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Fig. 5. The spectra of regions c and e. The lower spectra in (a) and (b) are the observed spectra downscaled by factors 
of 50 and 20, respectively. 



The mean value of the He mass fraction Y — 0.243 ± 0.008 
obtained for region c (see Table 0J, is consistent with the 
primordial 4 He mass fraction Y p = 0.244 ± 0.002, derived 
by extrapolating the Y vs O/H linear regression to O/H 
= (Izotov & Thuan lT§5g| >. and with Y p = 0.245 ± 0.002 
derived for the two most metal-deficient BCDs, known, I 
Zw 18 and SBS 0335-052 (Izotov et al. lT5§5)> . However, the 
He abundance in region c is likely slightly underestimated 
because of underlying stellar He I absorption. Due to the 
relatively small equivalent widths of the He I emission lines 
(Table El the effect of the underlying absorption may lead 
to an underestimate of the He + abundance by 1% - 2% 
in the case of the strongest He I A5876 line. This effect 
is larger for the weaker He I A4471 and A6678 emission 



lines. In region e only weak He I A3889 and A5876 emission 
lines are observed. Therefore, an accurate He abundance 
determination is not possible in that region. 

5. The stellar populations 

The low metallicity of HS 1442+4250 suggests that it may 
be a young nearby dwarf galaxy (Izotov & Thuan[l999). 
However, the morphological properties of the galaxy, such 
as the presence of multiple star-forming regions on the one 
hand, and a redder extended diffuse stellar component on 
the other hand, suggest that stars in HS 1442+4250 have 
been formed during several episodes of star formation. We 
study in this section the properties of the stellar popu- 
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Table 2. Observed (F(X)) and corrected (/(A)) fluxes and equivalent widths (EW) of emission lines in regions c and 
e. 



A (A) Ion 






region c 










region e 




F(A)/F(H/3) 


J(A)/I(H/3) 


EW(A) 


F(\)/F(UP) 


/(A)//(H/3) 


EW(k) 


3727 [O n] 


0.496 


±0.008 


0.535 


±0.009 


38.3 


±0.8 


1.477 


±0.089 


1.376 ±0.094 


24.4 ±1.8 


3750 H12 


0.020 


±0.004 


0.027 


±0.006 


1.6 


±0.4 










3771 Hll 


0.028 


±0.004 


0.035 


±0.006 


2.1 


±0.4 










3798 H10 


0.034 


±0.003 


0.042 


±0.006 


2.6 


±0.4 










3835 H9 


0.056 


±0.004 


0.065 


±0.006 


4.4 


±0.4 










3868 [Ne ill] 


0.405 


±0.007 


0.432 


±0.008 


32.1 


±0.7 


0.260 


±0.048 


0.242 ±0.048 


4.0 ±1.4 


3889 H8 + He i 


0.176 


±0.005 


0.193 


±0.007 


15.0 


±0.6 


0.070 


±0.033 


0.171 ±0.089 


1.2 ±1.0 


3968 [Nc in] + H7 


0.280 


±0.006 


0.302 


±0.007 


24.1 


±0.7 


0.144 


±0.034 


0.237 ±0.078 


2.5 ±1.0 


4026 He I 


0.015 


±0.004 


0.016 


±0.004 


2.0 


±0.5 










4101 US 


0.260 


±0.006 


0.277 


±0.007 


24.9 


±0.7 


0.174 


±0.032 


0.261 ±0.065 


3.1 ±1.0 


4340 H 7 


0.467 


±0.008 


0.486 


±0.008 


50.0 


±1.0 


0.374 


±0.039 


0.438 ±0.056 


7.4 ±1.4 


4363 [O in] 


0.134 


±0.004 


0.138 


±0.004 


14.2 


±0.7 


0.058 


±0.028 


0.055 ±0.028 


1.3 ±1.2 


4388 He i 


0.005 


±0.002 


0.005 


±0.002 


0.6 


±0.4 










4471 He i 


0.034 


±0.003 


0.035 


±0.003 


3.8 


±0.5 










4686 He n 


0.030 


±0.003 


0.030 


±0.003 


3.7 


±0.6 










4713 [Ar iv] + He I 


0.022 


±0.003 


0.022 


±0.003 


2.7 


±0.6 










4740 [Ar iv] 


0.014 


±0.002 


0.014 


±0.002 


1.8 


±0.5 










4861 H/3 


1.000 


±0.012 


1.000 


±0.012 


147.1 


±1.9 


1.000 


±0.063 


1.000 ±0.070 


25.7 ±2.5 


4922 He I 


0.010 


±0.003 


0.010 


±0.003 


1.4 


±0.6 










4959 [O in] 


1.769 


±0.019 


1.754 


±0.018 


250.0 


±2.6 


0.831 


±0.057 


0.774 ±0.057 


20.1 ±2.4 


5007 [O in] 


5.032 


±0.046 


4.976 


±0.046 


782.3 


±4.3 


2.489 


±0.130 


2.317 ±0.130 


62.5 ±6.9 


5876 He I 


0.099 


±0.003 


0.094 


±0.003 


22.7 


±1.3 


0.100 


±0.026 


0.093 ±0.026 


3.6 ±2.0 


6300 [O i] 


0.009 


±0.002 


0.008 


±0.002 


2.4 


±0.8 










6312 [S in] 


0.016 


±0.002 


0.014 


±0.002 


4.4 


±0.9 










6563 Ha 


3.012 


±0.028 


2.763 


±0.029 


881.8 


±6.8 


2.900 


±0.147 


2.738 ±0.162 


133.8 ±6.9 


6584 [N n] 


0.025 


±0.002 


0.023 


±0.002 


7.5 


±1.0 


0.070 


±0.019 


0.066 ±0.019 


3.3 ±2.0 


6678 He I 


0.028 


±0.002 


0.026 


±0.002 


8.6 


±1.1 










6717 [S n] 


0.046 


±0.002 


0.042 


±0.002 


13.9 


±1.2 


0.147 


±0.024 


0.137 ±0.025 


7.0 ±2.5 


6731 [S II] 


0.032 


±0.002 


0.029 


±0.002 


9.7 


±1.1 


0.095 


±0.020 


0.088 ±0.020 


4.6 ±2.1 


7065 He I 


0.023 


±0.002 


0.020 


±0.002 


7.8 


±1.1 










7136 [Ar m] 


0.043 


±0.002 


0.038 


±0.002 


15.8 


±1.4 










C(H/3) dex 






0.110±0.013 










0.000±0.063 




F(H/3) a 






3.85±0.03 










0.44±0.02 




EW(aba) A 






0.4±0.4 










1.45±1.0 





a in units 10 14 erg s 'cm 



lations in the extended LSB component and discuss the 
evolutionary status of HS 1442±4250. 

For this we use four different methods, three of them 
based on the spectroscopic data and the fourth one using 
photometric data. Two of these methods are based on the 
time evolution of equivalent widths of hydrogen Balmer 
line emission and absorption which are detected in al- 
most all regions along the slit. The other two methods are 
based on the comparison of observed with theoretical spec- 
tral energy distributions (SED) and broad-band colours. 
All the aforementioned methods are sensitive to the star 
formation history. The first two methods are extinction- 
insensitive, while the others depend on both the properties 
of the stellar population and the interstellar extinction. 
The use of all four methods allows us to put constraints 
on the star formation history and properties of the stellar 
population as well as to estimate the interstellar extinc- 



tion. These methods are described in detail in Guseva et 

ai. (EnnnEonsii- 

The observed (V — I) colour of the brightest region c is 
very blue ~ -0.5 (Fig.[5J;), and cannot be reproduced by a 
stellar population of any age. This is because of the large 
contribution of very strong oxygen and hydrogen emission 
lines, as evidenced by their large equivalent widths (Table 

and gaseous continuum. Gaseous emission is also an 
important contributor to the total light of regions I and 
II. However, in other regions the equivalent widths of the 
emission lines are small (Table . Therefore, for those 
regions we will not take into account ionized gas emission. 

5.1. Age from Balmer nebular emission lines 

The equivalent widths of hydrogen nebular emission lines 
are usually used as an age indicator of star-forming re- 
gions, in which O and early B stars are still present. The 
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equivalent widths of these lines for an instantaneous burst 
decrease sharply after ~ 10 Myr (e.g., Schaerer & Vacca 
1998). However, in the case of continuous star formation 
this method can be succesfully used for age determination 
in the 0-10 Gyr range (Guseva et al. [2001). 

For the age determination we use the two brightest hy- 
drogen emission lines, Ha and H/3. Their fluxes and equiv- 
alent widths, measured in the spectra of the extracted LSB 
regions are shown with errors in Table 

Because the H/3 emission line is narrower than the ab- 
sorption line in these regions and does not fill the absorp- 
tion component, its flux was measured using the contin- 
uum level at the bottom of the absorption line. This level 
has been chosen by visually interpolating from the absorp- 
tion line wings to the center of the line. 

The extinction coefficients C(H/3) in regions with de- 
tectable Ha and H/3 emission lines are derived from the 
observed Ha/H/3 flux ratios. The theoretical recombina- 
tion Ha/H/3 flux ratio of 2.8 is adopted, which is typical 
for low-metallicity H II regions. No correction for the ab- 
sorption line equivalent widths has been made. The ex- 
tinction coefficients C(H/3) are shown in Table 21 Note, 
that C(H/3) for region / is uncertain because of the low 
equivalent width of the H/3 emission line which is compa- 
rable to the equivalent width of the absorption line. 

The dependence of the Ha and H/3 emission line equiv- 
alent widths on age for a heavy element mass fraction 
Zq/20 has been calculated for two limiting star formation 
histories, the case of an instantaneous burst, and that of a 
continuous constant star formation. The calculations are 
done with the galactic evolution code PEGASE.2 (Fioc 
& Rocca- Volmerange 11997|) . The time evolution of Ha in 
the case of an instantaneous burst is shown in Fig. 01 by 
a thick solid line. We next consider the case of continuous 
star formation with a constant star formation rate (SFR). 
The temporal evolution of the equivalent widths of the 
H/3 and Ha emission lines is calculated using the model 
equivalent widths of hydrogen emission lines and SEDs for 
an instantaneous burst (Fioc & Rocca- Volmerange[[997). 
The temporal dependences of the Ha emission line equiv- 
alent width are shown in Fig. by different thin lines 
for continuous star formation starting at time t- u as de- 
fined by the abscissa value, and stopping at tf (the curves 
in Fig. are labeled by log tf). Time is zero now and 
increases to the past. 

To compare theoretical predictions with observed data 
we consider in the following regions a and / only, as rep- 
resentative of the whole LSB component. The positions of 
the measured EW(Ra) on the model curve for continuous 
star formation with tf = 8 Myr are shown in Fig. [HJi by a 
triangle for region a and by a star for region /. In this case 
the data are consistent with star formation starting not 
earlier than ti ~ 20 Myr ago for region a and ~ 200 Myr 
ago for region /. If instead continuous star formation stop- 
ping at tf = is considered, then the observed EW(R/3) 
and EW(Ha) are consistent with models in which star for- 
mation starts at time ti ~ 1 — 4 Gyr ago for region a and 
ti ^> 10 Gyr ago for region /. Hence, we conclude that for 



region /, models in which stars are continuously forming 
until now (tf — 0) with a constant star formation rate are 
inconsistent with the observed EW(Ha) and EW(HP). 

5.2. Age from Balmer stellar absorption lines 

Hydrogen Balmer absorption lines are detected along the 
slit in a large part of HS 1442+4250. This allows us to 
estimate the age of stellar populations by comparing the 
observed hydrogen absorption line equivalent widths with 
theoretical ones. Gonzalez Delgado et al. (1999J have cal- 
culated the temporal dependence of the Balmer absorp- 
tion lines considering both instantaneous and continuous 
star formation. Their models for an instantaneous burst 
predict a steady increase of the equivalent widths with 
time for ages ranging from 1 Myr to 1 Gyr, reaching max- 
imum values EW (B.-y) ~ 12A and EW(RS) ~ 15A for a 
metallicity Z — Zq/20 (Fig. Et>) ■ For solar metallicity, the 
predicted maximum EWs are larger by ~ lA. However, 
Gonzalez Delgado et al. (1999) do not extend calculations 
to ages > 1 Gyr, when the equivalent widths of the ab- 
sorption lines decrease with age (Bica & Alloin^986a). 

We compared the model predictions of the Hi5 absorp- 
tion line equivalent width for an instantaneous burst by 
Gonzalez Delgado et al. (1999) (Fig. 9 in Guseva et al. 
2003) with the empirical calibration by Bica & Alloin 
(1986a) (Fig. |HJd) and with available observational data 
for a large sample of different objects: open and globu- 
lar stellar clusters (Bica & Alloin [T"986alll986b|) . nuclei of 
normal elliptical and spiral galaxies (Bica & Alloin[1987 
Bica QUO Schmidt et 81. 115551115551 Saraiva et al.HEHJ) 
with different ages and metallicities. No stellar clusters 
and galaxies with £W(Hy) > 9.6A and EW(HS)> 12A 
were found in the samples considered. Hence, the models 
by Gonzalez Delgado et al. I|1999|l at ages ~ 1 Gyr appar- 
ently overestimate the equivalent widths of the absorption 
lines. 

Therefore we only use an empirical calibration of 
Balmer absorption line equivalent widths versus age by 
Bica & Alloin (1986a) for ages ranging between 1 Myr 
and 16.5 Gyr. This calibration is based on integrated spec- 
tra of 63 star clusters with known ages, metallicities and 
reddenings. Later, a calibration of the H7 and H<5 absorp- 
tion line equivalent widths was derived by Schmidt et al. 
(1995) to apply to starburst events in dwarf galaxies. The 
dependence of the equivalent width of the H5 absorption 
line on age for an instantaneous burst from the data of 
Bica & Alloin (1986a) is shown in Fig. 03 by thick solid 
line. 

We also consider the case of continuous star formation. 
For this we assume that stars are forming with a constant 
star formation rate starting at time t\ and stopping at tf. 
We use the equivalent widths of hydrogen absorption lines 
from the empirical calibration by Bica & Alloin (1986a) 
and SEDs by Fioc & Rocca- Volmerange : 1 9iT7l) for instan- 
taneous bursts to calculate the temporal evolution of the 
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Table 4. Fluxes, equivalent widths of Ha and H/3 emission lines and the extinction coefficient C(H/3) in LSB regions. 



Region 


Distance 3 


Aperture b 


F(Hq) c 




F(H/3) C 


EW(Hf3) d 


C(H0) 


I 


-24.4 


2.0x4.1 


8.0 ±0.5 


1200.0 ±2.5 


3.0 ±0.5 


32.8 ±0.8 


0.0 


a 


-12.4 


2.0x7.6 


50.4 ±0.6 


78.9 ±1.8 


17.9 ±0.5 


14.6 ±0.7 


0.00 ±0.07 


b 


-5.8 


2.0x4.1 


22.0 ±0.6 


36.5 ±0.8 


7.2 ±0.6 


7.3 ±0.5 


0.12 ±0.07 


d 


8.0 


2.0x5.5 


53.9 ±0.3 


52.4 ±2.5 


17.1 ±0.6 


9.9 ±0.5 


0.15 ±0.07 


f 


23.9 


2.0x4.1 


13.6 ±0.6 


16.8 ±1.8 


2.8:±0.5 


3.1 ±0.7 


0.31:±0.07 


fJ 


28.4 


2.0x4.8 


17.4 ±0.6 


31.7 ±0.8 








h 


34.2 


2.0x6.9 


5.5 ±0.3 


12.4 ±1.9 








II 


44.2 


2.0x6.2 


13.2 ±0.3 


634.4 ±2.5 


4.9 ±0.6 


32.4 ±0.5 


0.0 



a distance from region c in arcsec. Negative and positive values correspond to regions located respectively to the southwest and 
northeast from region c. 

b aperture x x y where x is the slit width and y the size along the slit in arcsec. 
c in units 10 -16 erg s -1 cm -2 . 
d in A. 



Table 5. Equivalent widths of H7 and H<5 absorption lines in LSB regions. 



Region 


Distance 3 


Aperture b 




£W(H 7 ) c ' d 


EW(my^ A 




a 


-12.4 


2.0x7.6 


9.0:±0.6 




9.8:±0.6 




b 


-5.8 


2.0x4.1 


8.6 ±0.4 


6.6 ±0.5 


8.9 ±0.4 


6.7±0.4 


d 


8.0 


2.0x5.5 


7.3 ±0.5 


5.3 ±0.5 


7.8 ±0.5 


5.0±0.5 


f 


23.9 


2.0x4.1 


6.2 ±0.7 


5.1 ±0.9 


6.6 ±0.7 


5.1±0.7 


9 


28.4 


2.0x4.8 


7.2 ±0.4 


4.8 ±0.9 


7.8 ±0.4 


5.1±0.4 


h 


34.2 


2.0x6.9 


7.4 ±0.5 


5.3 ±0.5 


7.6 ±0.5 


7.4±0.5 



a distance from region c in arcsec. Negative and positive values correspond to regions located respectively to southwest and 
northeast from region c. 

b aperture x x y where x is the slit width and y the size along the slit in arcsec. 
c in A. 

Equivalent widths are measured using Gaussian fitting of the lines. 

Equivalent widths are obtained within the same wavelength intervals Ao = 4318 - 4364A and 4082 - 4124A, respectively, for 
H7 and H<5 as those used by Bica & Alloin (1986a ). 



equivalent widths of hydrogen absorption lines for contin- 
uous star formation. 

The results are given in Fig. 03. The temporal depen- 
dences of the equivalent width of the US absorption line 
are shown by different thin lines for continuous star forma- 
tion starting at t\, as defined by the abscissa value, and 
stopping at U (the curves in Fig. 03 are labeled by log 
tf). The equivalent width of the H<5 absorption line in the 
spectrum of the stellar population formed between t\ and 
tf corresponds in Fig. 03 to the EW at time t\. 

Stellar hydrogen absorption lines have been detected 
in all LSB regions. In all regions except for region a, we 
measure the equivalent widths of the US and H7 absorp- 
tion lines. The H7 absorption line in region a has not been 
used because of the strong contamination by nebular emis- 
sion. The contribution of the nebular H7 and US emission 
lines to the corresponding absorption lines in the spectra 
of regions a, 6, d, /, g and h has been removed using the 
IRAF software package. For the northeastern regions /, 
g, and h, the contamination of the US and H7 absorption 
lines by nebular emission is small, resulting in a correction 
which is within the errors of the EWs. 

Table Q3 lists the equivalent widths with errors of the 
H7 and US absorption lines. The equivalent widths shown 



in cols. (4) and (5) are measured using Gaussian fitting 
of the lines, while the EWs in cols. (6) and (7) marked 
with the subscript "B&A" are obtained within the same 
wavelength intervals A = 4318 - 4364A and 4082 - 4124A, 
respectively, for H7 and H<5 as those used by Bica & Alloin 

The positions of the observed EW(HS) on the model 
curve for continuous star formation with tf = 8 Myr are 
shown in Fig. ^jp by triangles for region a and by a star for 
region /. The observed EW(H5) and EW(Hj) are consis- 
tent with those of stellar populations forming continuously 
from ~ 40-250 Myr ago to 8 Myr ago if the model calcu- 
lations by Gonzalez Delgado et al. l|1999f) are used. These 
models are shown in Fig. 9 of Guseva et al. (2003). In the 
case of the empirical calibration by Bica & Alloin (1986a) 
the observed EW{US) (Fig. \flp) and EW(Rj) are best 
reproduced by a stellar population forming continuously 
with a constant star formation rate starting at i; ~ 50 
Myr for region / and ~ 1 Gyr for region a and stopping at 
tf = 8 Myr. 

Comparison of Figs. HO* and[U> shows that the ages of 
the oldest stars contributing to the light of region / (stars) 
derived from the emission and absorption line equivalent 
widths are consistent and lie in the narrow interval be- 
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log tj (yr) log t l (yr) 

Fig. 6. (a) Temporal evolution of the equivalent width of the Ha emission line for an instantaneous burst (thick solid 
line) and a heavy element mass fraction Z — Zq/20 calculated with the PEGASE.2 code (Fioc & Rocca-Volmerange 
1997). Other different lines show model predictions for continuous star formation with constant SFR, starting at the 
time U ago, defined by the abscissa, and stopping at tf, with tf — (thin solid line) and tf changing from 5 Myr to 30 
Myr ago. Each line is labeled by the respective log tf (tf in yr). The positions of the observed EW(Ha) on the modeled 
curve are shown for region a by triangle and for region / by star for the case of continuous star formation with if = 
8 Myr. (b) Equivalent width of the H<5 absorption line vs age for stellar clusters (Bica & Alloin 1986a) (thick solid 
line). The predicted EW(HS) are also shown for continuous star formation with constant SFR starting at time i; ago, 
defined by the abscissa, and stopping at if, with tf — (thin solid line) and tf from 5 Myr to 6.3 Gyr. Each line is 
labeled by the respective log tf (tf in yr). The positions of the observed EW(RS) on the modeled curve for region a are 
shown by triangles and the position of the observed EW(HS) for region / is shown by star for the case of continuous 
star formation with tf = 8 Myr. The equivalent widths have been measured in the same wavelength intervals as those 
used by Bica & Alloin l|1986a|) . 



tween 50 - 150 Myr. This implies that continuous star 
formation stopping at tf = 8 Myr is a reasonable scenario 
for region /. However, the ages of the stellar population in 
region a derived from Figs.^ andlHJ) (triangles) are not in 
agreement. The age of the oldest stars contributing to the 
light of this region as derived from the emission lines is ~ 
25 Myr, while the age derived from the absorption lines us- 
ing the Gonzalez Delgado et al. ■ ] 1 99 9 ■ and Bica & Alloin 
( 1986a) calibrations is ~ 300 Myr and ~1 - 2 Gyr, respec- 
tively. Evidently, continuous star formation stopping at tf 
= 8 Myr does not reproduce the observed properties of re- 
gion a. A better agreement is achieved if star formation in 
this region continues until now (i.e. tf = 0, thin solid lines 
in Fig. [(J- andEJ}). Then, within the errors, the observed 
data gives ages ~500 Myr (using Gonzalez Delgado et al., 
1999) and ~1 Gyr (Fig.Eb) for the oldest stars, which is in 
better agreement with the age estimation ^2 - 3 Gyr ob- 
tained from the equivalent widths of emission lines. Note, 
that choosing tf = is a limiting case which gives a max- 
imum age of the old stellar population in a scenario with 
a constant star formation rate. 



5.3. Age from the spectral energy distribution 

We next consider the age of the stellar population derived 
from the comparison of the observed and theoretical SEDs. 
The shape of the observed spectrum reflects not only the 
properties of the stellar population but also reddening ef- 
fects. Therefore only a combination of the SED method 
with those discussed in Sect. 15.11 and 15.21 allows us to si- 
multaneously estimate age and interstellar extinction for 
the LSB regions in HS 1442+4250. 

To fit the observed SEDs, we use the galactic evolu- 
tion code PEGASE.2 (Fioc & Rocca-Volmerange lTMTjl to 
produce a grid of theoretical SEDs for an instantaneous 
burst of star formation with ages ranging between and 
10 Gyr, and a heavy element mass fraction Z = Zq/20. 
An initial mass function with a Salpeter slope (a = -2.35), 
and upper and lower mass limits of 120 Mq and 0.1 Mq 
are adopted. The contribution of gaseous emission to the 
total emission of regions a, b, d, /, g and h is small and has 
not been taken into account in the SED calculations. 

A reliable determination of reddening can only be 
made for the two brightest H II regions c and e because 
of their many strong hydrogen emission lines (Table0). In 
other regions, only Ha and Hf3 emission lines are present. 
The latter line is very weak in regions /, g and h. The 
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extinction coefficient C(H/3) obtained from the Balmcr 
decrement in the different regions ranges from to 0.15. 
C(H/3) is 0.11 in the brightest region c. 

We consider only the properties of region /, taken to 
be representative of the LSB region. We showed in Sect. 
I5.1l and l5.2l that the observed equivalent widths of the hy- 
drogen emission and absorption lines in this region can be 
reproduced by a stellar population formed continuously 
with a constant SFR, starting at t\ ~ 100 Myr and stop- 
ping at tf = 8 Myr. However, the observed SED of region 
/ cannot be reproduced by a synthetic SED with such a 
young stellar population (lower spectra in FigEJi) without 
assuming a non- negligible reddening. In particular, the ob- 
served SED can be fitted by a synthetic SED of a stellar 
population formed continuously between 10 and 50 Myr 
with a constant SFR, if there is an extinction C(H/3) = 
0.29 (upper spectra in Fig. EJi). 

Alternatively, the observed SED of region / at A <; 
3900A can also be well fitted by a synthetic SED of a 
stellar population formed continuously between and 10 
Gyr, if the SFR during the last 50 Myr was enhanced by 
a factor b of 10 and no reddening is present (SED labeled 
"b=10" in Fig. |7Jd). However, the fit is not satisfactory 
at shorter wavelengths. Additionally, this model does not 
reproduce the observed EW (Ha) and EW(H(3). 

We also consider models with star formation occurring 
in two episodes separated by a quiescent period. A model 
with star formation occurring with a constant SFR during 
(1) 1-10 Gyr and (2) 10 - 50 Myr with comparable 
masses of stellar populations formed in each period (not 
shown in Fig. 0) is able to reproduce the observed EWs 
of the hydrogen lines (EW(Ua) = 111, EW{K0) = 2A, 
£W(H 7 ) = -5. 7 A and EWQ1S) = -6.7A). However, in 
that case the observed SED is again not well fitted at A 
< 3900A. All the observational data for region / are best 
fitted by a model in which the stellar population formed 
continuously with a constant SFR in two periods 1) 0.5 - 
2 Gyr and 2) 8-30 Myr ago and no reddening is present 
(Fig. Ek). If some extinction is present then ages will be 
smaller. 

In summary, we conclude that only young and 
intermediate-age stellar populations can reproduce the ob- 
served properties of region /. 

5.4. Age from the colour distribution 

We derived V and I surface brightness and colour distri- 
butions for the regions covered by the spectroscopic ob- 
servations and compared them with predictions from our 
population synthesis modeling. The results of this com- 
parison are shown in Fig. The predicted colours ob- 
tained from convolving the theoretical SEDs with the ap- 
propriate filter bandpasses are shown by different symbols. 
The transmission curves for the Johnson V and Cousins I 
bands are taken from Bessell (1990). The zero points are 
from Bessell et al. (JTHBJ) • 



The contribution of the gaseous emission to the total 
brightness is small in the LSB component. Therefore, we 
do not take it into account. The colours of regions a and 
b are well fitted by those of a stellar population formed 
continuously with a constant SFR in two periods of star 
formation: (1) 100 - 400 Myr and (2) 4-15 Myr ago. The 
colours of the northeastern part of the galaxy excluding 
the brightest region c are best fitted by those of a stellar 
population forming continuously with a constant SFR in 
two periods: (1) 0.5 - 2 Gyr and (2) 8-30 Myr ago. In 
all cases we adopt C(H/3)=0, which gives upper limits to 
the ages. 

Since the contribution of ionized gas emission to the 
total light of the brightest H n region (region c) is high, the 
theoretical SED for this region has been constructed using 
a 4 Myr stellar population SED for a heavy element mass 
fraction Z = Zq/20, and adding the gaseous continuum 
SED and observed emission lines (see Guseva et al. 2001 
for details). The predicted (V — I) colour of region c is 
reddened adopting C(H/3) = 0.11 (Table |2J. The observed 
colour (V — I) ~ -0.5 mag of this region is bluer than 
that of a 4 Myr stellar population (labeled by an asterisk 
in Fig. |5p) implying a large contribution of ionized gas 
emission (labeled by a triangle in Fig. Eb)- 

In Fig. 13; we show with filled circles the modeled 
colours of all regions in the LSB component of the galaxy 
and the total colour of region c. The agreement between 
the observed and synthetic (V — I) colours is very good. 
The reddest colours (V — I) ~ 0.63 mag in HS 1442+4250 
are found in the outer parts of the galaxy at \iy > 24 mag 
arcsec -2 (Fig. |3}. No spectroscopic data are available for 
these regions. Therefore, we can only use photometric data 
to estimate the age of the stellar population in those re- 
gions. Again, the maximum age depends on the adopted 
star formation history. If an instantaneous burst model is 
assumed then the age of the population with (V— I) = 0.63 
mag is ~1 Gyr. In the case of continuous star formation 
the same colour can be explained by a stellar population 
formed with a constant star formation rate between and 
5 Gyr. We note, however, that the age estimates for the 
outer part of HS 1442+4250 are uncertain due to large 
errors in the (V — I) colour and likely contamination of 
the LSB light by red background/foreground objects. 

6. Conclusions 

We present a detailed photometric and spectroscopic 
study of the low-metallicity dwarf irregular galaxy HS 
1442+4250. Broad-band V and / images and spectra in 
the optical range have been obtained with the 2.1m and 
4m Kitt Peak telescopes respectively. The main conclu- 
sions of this study can be summarized as follows: 

1. HS 1442+4250 is a low-metallicity nearby (D 

I '2.1 Mpc) dwarf irregular galaxy with a chain of H 

II regions arranged along an elongated low-surface- 
brightness (LSB) component. The LSB component is 
well fitted by an exponential profile with a scale length 
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a > 270 pc in its outer part, flattening for small radii 
(R* <, 2a). The observed (V - I) colour of the bright- 
est H II region is very blue ~ -0.5 mag due to the 
combined effect of a young stellar population and ion- 
ized gas emission. The colour of the LSB component 
is constant and equal to ~ 0.4 mag in a large part of 
the galaxy. 

2. The oxygen abundance for the brightest H n region 
c is 12 + log(0/H) = 7.63 (Z = Z G /19). The neon- 
to-oxygen abundance ratio for this region, log Ne/O 
= -0.75, is in good agreement with the mean ratio 
derived from previous studies of dwarf galaxies (e.g., 
Izotov & Thuan ■ The nitrogen-to-oxygen abun- 
dance ratio log N/O = -1.44 is ~ 0.15 dex larger than 
the N/O ratios obtained by Thuan et al. {1.995) and 
Izotov & Thuan l|1999J) for the most metal-deficient 
blue compact dwarf (BCD) galaxies with Z <, Z Q /20. 

3. The 4 He mass fraction Y = 0.243 ± 0.008 derived for 
the brightest H n region is in good agreement with 
the primordial 4 He mass fraction Y p = 0.244 - 0.245 
derived by Izotov & Thuan (1998) and Izotov et al. 

{B2EJ. 

4. We use four methods to estimate the age of the stellar 
population in the LSB component of HS 1442+4250. 
Different histories of star formation are considered. 
The spectroscopic data such as the equivalent widths 
of the hydrogen Ha and H/3 emission lines, and of the 
hydrogen H7 and US absorption lines, the spectral en- 
ergy distributions and the photometric colours are re- 
produced quite well by models with only young and 
intermediate-age stellar populations. It is unlikely that 
a significant old stellar population contributes much 
to the observed light in the inner regions. However, 
in the outermost regions of the LSB component, with 
(V — I) ^0.63, the presence of older populations is not 
excluded. For these regions only photometric data are 
available. The observed colour can be explained by a 
stellar population continuously forming between and 
5 Gyr with a constant SFR. However, because of large 
photometric errors and possible confusion with back- 
ground sources, those age estimates are uncertain. 

Acknowledgements. N.G.G. and Y.I.I, acknowledge DFG grant 
436 UKR 17/2/02 and Y.I.I, is grateful for the Gauss profes- 
sorship of the Gottingen Academy of Sciences. They are also 
grateful for Swiss SCOPE 7UKPJ62178 grant and for hospi- 
tality at Gottingen Observatory. Y.I.I, and T.X.T. acknowl- 
edge partial financial support through NSF grant AST-02- 
05785. Research by P.P. and K.J.F. has been supported by 
the Deutsches Zentrum fur Luft- und Raumfahrt e.V. (DLR) 
under grant 50 OR 9907 7. K.G.N, thanks support from the 
Deutsche Forschungsgemeinschaft (DFG) grants FR 325/50-1 
and FR 325/50-2. 

References 

Anders, E., & Grevesse, N. 1989, Geochim.Cosmochim.Acta, 
53, 197 

Bessell, M. S. 1990, PASP, 102, 1181 



Bessell, M. S., Castelli, F., & Plez, B. 1998, A&A, 333, 231 
Bica, E. 1988, A&A, 195, 76 
Bica, E., & Alloin, D. 1986a, A&AS, 66, 171 
Bica, E., & Alloin, D. 1986b, A&A, 162, 21 
Bica, E., & Alloin, D. 1987, A&AS, 70, 281 
Binggeli, B., & Cameron, L. M. 1991, A&A, 252, 27 
Brocklehurst, M. 1971, MNRAS, 153, 471 
Caon, N., Capaccioli, M., & D'Onofrio, M. 1993, MNRAS, 265, 
1013 

Fioc, M., & Rocca-Volmerange, B. 1997, A&A, 326, 950 
Fricke, K. J., Izotov, Y. I., Papaderos, P., Guseva, N. G., & 

Thuan, T. X. 2001, AJ, 121, 169 
Garnett, D. R. 1992, AJ, 103, 1330 

Gonzalez Delgado, R. M., Leitherer, C, & Heckman, T. M. 

1999, ApJS, 125, 489 
Guseva, N. G., Izotov, Y. I., & Thuan, T. X. 2000, ApJ, 531, 

776 

Guseva, N. G., Izotov, Y. I., Papaderos, P., et al. 2001, A&A, 
378, 756 

Guseva, N. G., Papaderos, P., Izotov, Y. I., et al. 2003, A&A, 
in press 

Izotov, Y. I., & Thuan, T. X. 1998, ApJ, 500, 188 
Izotov, Y. I., & Thuan, T. X. 1999, ApJ, 511, 639 
Izotov, Y. I., Thuan, T. X., & Lipovetsky, V. A. 1994, ApJ, 
435, 647 

Izotov, Y. I., Lipovetsky, V. A., Chaffee, F. H., et al. 1997, 
ApJ, 476, 698 

Izotov, Y. I., Chaffee, F. H., Foltz, C. B., et al. 1999, ApJ, 527, 
757 

Kniazev, A. Y., Pustilnik, S. A., & Ugryumov, A. V. 1998, 

Bulletin Spec. Astrophys. Obs., 46, 23 
Kraan-Korteweg, R. C. 1986, A&AS, 66, 255 
Landolt, A. U. 1992, AJ, 104, 340 
Makarova, L. 1999, A&A, 139, 491 

Papaderos, P., Loose, H.-tL, Thuan, T. X., & Fricke, K. J. 

1996a, A&AS, 120, 207 
Papaderos, P., Loose, H.-H., Fricke, K. J., & Thuan, T. X. 

1996b, A&A 314, 59 
Popescu, C. C, & Hopp, U. 2000, A&AS, 142, 247 
Popescu, C. C, Hopp, U., Hagen, H. J., & Elsasser, H. 1996, 

A&AS, 116, 43 

Pustilnik, S. A., Engels, D., Ugryumov, A. V., et al. 1999, 

A&AS, 137, 299 
Sanduleak, N., & Pesch, P. 1982, ApJ, 258, Lll 
Saraiva, M. F., Bica, E., Pastoriza, M. O, & Bonatto, C. 2001, 

A&A, 376, 43 
Schaerer, D., & Vacca, W. D. 1998, ApJ, 497, 618 
Schmidt, A. A., Alloin, D., & Bica, E. 1995, MNRAS, 273, 945 
Schmidt, A. A., Bica, E., & Dottori, H. A. 1989, MNRAS, 238, 

925 

Sersic, J. L. 1968, Atlas de Galaxias Australes (Cordoba: Obs. 
Astron.) 

Stasihska, G. 1990, A&AS, 83, 501 

Thuan, T. X., Izotov, Y. I., & Lipovetsky, V. A. 1995, ApJ, 
445, 108 

Thuan, T. X., Izotov, Y. I., & Foltz, C. B. 1999, ApJ, 525, 105 
Tifft, W. O, Kirshner, R. P., Gregory, S. A., & Moody, J. W. 

1986, ApJ, 310, 75 
Vennik, J., Hopp, U., Kovachev, B., et al. 1996, A&AS, 117, 

216 

Vennik, J, Hopp, U, & Popescu, C. C. 2000, A&AS, 142, 399 



14 



N.G. Guseva et al. 



6.0 



4.0 



region / 



lC(H/3)=0.29 



continuous 8F 
10 - 50 Myr 



2.0 



EW =3 _ 
EW H =-5.5 
EW„,,= -6.7. 




2.0 



C(H0) = O 



two periods of SF 



(c) 



1) 500Myr-3Gyr 

2) 8 - 30 Myr 





5 10 15 30 25 



R (arcsec) 

Fig. 4. (a) Surface brightness profiles (SBPs) of HS 
1442+4250 in V (filled circles) and I (open circles) de- 
rived using the method iii (Papaderos et al. I1996afl . The 
modeled surface brightness distribution of the LSB com- 
ponent in V assuming a modified exponential distribution 
by Papaderos et al. H199fia|) is shown by the thick-grey 
curve. This model implies a central surface brightness for 
the LSB component ss 1.3 mag fainter than the value pre- 
dicted by extrapolation of the outer exponential slope of 
the SBPs to R* = 0". The surface brightness distribution 
of the light in excess of the LSB component (crosses) is 
due to the emission of star-forming regions arranged along 
the major axis of the galaxy (starburst component) (see 
Figs.HandEH. (b) (V-I) colour profile of HS 1442+4250 
computed by subtraction of the SBPs displayed in the up- 
per panel. The thick-grey curve shows the (V — I) colour 
distribution obtained by subtracting the adopted models 
of the LSB component intensity distribution in V and I. 
Both profiles suggest a slow colour increase of 0.3 mag 
kpc -1 with radius and a mean (V — I) colour of 0.63 ± 
0.05 mag for R* > 15". The effective radius r c g of the V 
SBP and the isophotal radii P25 and E25 of the starburst 
and LSB components at 25 V mag arcsec -2 are indicated. 
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Fig. 7. Spectrum of region / (a) — (c) on which synthetic 
SEDs are superposed. They are calculated for stellar pop- 
ulations formed continuously within different time inter- 
vals. The predicted equivalent widths of hydrogen emis- 
sion and absorption lines are shown for each model, (a) 
Synthetic SED corresponds to a stellar population formed 
continuously between 10 and 50 Myr ago with a constant 
SFR. The lower SED is superposed on the spectrum un- 
corrected for extinction and the upper SED is superposed 
on the spectrum corrected for extinction with C(H/3) = 
0.29. (b) Synthetic SED labeled b = SFR(i < t b )/SFR(f 
> th) = 1 with <b = 50 Myr, corresponding to a stel- 
lar population formed continuously between and 10 Gyr 
with a constant SFR. The synthetic spectrum labeled b 
= 10 corresponds to a stellar population formed contin- 
uously between and 10 Gyr with a SFR enhanced by 
10 times during t < t^. (c) Synthetic SED, calculated for 
star formation with a constant SFR occurring in two pe- 
riods: 1) 500 Myr - 2 Gyr ago and 2) 8 - 30 Myr ago. The 
extinction coefficient is C(H/3) = for (b) and (c). 



